Multi-stable regime and intermediate solutions in a nonlinear saturable coupler 
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We show, theoretically and experimentally, the existence of a multi-stable regime in a nonlinear 
saturable coupler. In spite of its simplicity, we found that this model shows generic and fundamental 
properties of extended saturable lattices. The study of this basic unit becomes crucial to under- 
stand localization mechanisms and dynamical properties of extended discrete nonlinear saturable 
systems. We theoretically predict the regions of existence of intermediate solutions, and experi- 
mentally confirm it by observing a multi-stable propagation regime in a LiNbC>3 saturable coupler. 
This constitutes the first experimental evidence of the existence of these unstable symmetry-broken 
stationary solutions. 
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I. INTRODUCTION 

Nonlinear discrete systems appear in several branches 
of science and have found a fruitful field of development 
and realistic implementation during the last years [H-[3|- 
Many results obtained in very different physical settings 
can be extrapolated to other areas of research generating 
a broader and deeper scientific impact. Different tech- 
niques and methods to study such systems experimen- 
tally, for example in the context of photonic lattices, have 
been developed with many possibilities to change and 
control the key parameters. Here, merely Kerr-like (cu- 
bic) nonlinear systems have been a main subject of the- 
oretical and experimental research. As a result, the cor- 
roboration of several former theoretical predictions and 
different new findings have been performed. However, a 
different type of nonlinearity has opened new challenges 
with new interesting dynamical properties, the so-called 
"saturable nonlinearity" . On one hand, from a dynam- 
ical point of view, this type of nonlinearity allows for a 
more complex and richer phenomenology than typical cu- 
bic systems For example, in a saturable nonlinear 
regime an exchange of stability properties between funda- 
mental solutions is allowed, promoting improved mobil- 
ity for high power solitons. This is certainly opposite to 
the phenomenology observed for cubic systems [13, [H| , 
where an increment in power induces localization only. 
On the other hand, only very little experimental realiza- 
tions showing the specific behavior of saturable systems 
have been performed. So far, mainly gradual changes 
from the behavior of Kerr-like systems were observed; 
for example, suppression of modulation instability [Hj], 
stabilization of discrete vector solitons [l3j], and higher- 
order gap solitons [lH ]. 

A nonlinear coupler (dimer) is the simplest discrete 
system where there are just two identical nonlinear 
waveguides (sites) - placed in close proximity - which 
evanescently interact. Thirty years ago, Jensen [TH ] 
showed the main features for the corresponding cubic 



system; i.e., a periodic exchange of light for small power 
and high transmission (localization) for larger powers. In 
Ref. [16( | authors explored the dynamics and stationary 
behavior of a cubic coupler which presents only one bi- 
furcation point for stationary solutions and no exchange 
of stability properties [similar to larger one-dimensional 
(ID) cubic systems]. On the other hand, Ref. [l7| shows 
the appearance of an extra bifurcation point, change of 
stability properties and richer dynamics when consider- 
ing a saturable nonlinearity. However, only very recently, 
the concept of "intermediate solutions" (IS) was intro- 
duced to explain theproperties of extended saturable ID 
and 2D lattices LD-lal (this concept was introduced be- 
fore for other models [ID, HH). These kind of "unstable" 
symmetry-broken solutions appear when two fundamen- 
tal modes are simultaneously stable (in other settings, 
the IS can also be stable [191 - 4 211 ) ) . Therefore, it becomes 
natural to formulate the question about the minimum 
number of sites - in a saturable array - for which this 
phenomenology emerges. In the present work, we will 
show that the fundamental saturable properties are al- 
ready present in a system of just two waveguides. More- 
over, for the first time to our knowledge, we show an ev- 
idence for the existence of IS by observing a multi-stable 
propagation in an experiment performed in iron doped 
LiNb03 samples. 



II. MODEL 

The propagation of light in a system composed of two 
identical weakly coupled waveguides, with a defocusing 
saturable nonlinearity, can be described as follows 
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where u n represents the light amplitude at site n , 7 = 
7/V > corresponds to the strength of the defocusing 
nonlinearity (7) with respect to the coupling coefficient 
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(V) between the two sites, and z describes the normalized 
propagation distance along the waveguides. Model ([T]) 
possesses two conserved quantities, the Power 



p= M 2 + M 2 



(2) 



and the Hamiltonian 



H=(u 2 u* 1 +u* 2 u l )+ 1 \n[(l + \u l \ 2 ){l + \u 2 \ 2 )} . (3) 

Stationary solutions of model ([TJ) have the form u n (z) — 
M„exp(iAz), where A represents the spatial frequency. 
First of all, we look for linear solutions of model ([T]) that, 
due to the saturation, exist in two different regions 0-Q ■ 
In a low power regime, we find that there are two solu- 
tions with frequencies "7 + 1" and "7 — 1" . For higher 
powers, the nonlinear response vanishes and frequencies 
become and "—1". The spatial profiles for these 

two modes are equal in both regimes: the symmetric 
(u\ = u 2 ) and the antisymmetric (ui = —u 2 ) modes, 
respectively. 
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Figure 1: (Color Online) (a) P vs A diagram for 7 = 10. 
The symmetric, stable and unstable antisymmetric, asym- 
metric, and intermediate solutions are plotted in black, full 
and dashed blue lines, red, and red-dashed line, respectively. 
Vertical lines indicate linear frequencies, (b) Antisymmetric, 
asymmetric, and intermediate mode profiles, for P = 20. 



III. NONLINEAR MODES 

Now, we look for general nonlinear solutions of the 
form: ui = A and u 2 — aA, where A is a positive am- 
plitude and a describes the ratio between these two site 
amplitudes. We found that the symmetric (a = 1) and 
the antisymmetric (a — —1) modes are also solutions in 
the nonlinear regime. They bifurcate from the low power 
linear solutions and diverge when approaching the high 
power linear modes: 



P — 2 

J sym ^ 



( 7 + l)-A 



A- 1 



and P ant = 2 



( 7 -l)-A 



A + l 



In Fig. Ufa) these two families are plotted with thick- 
black and blue lines, respectively, including their corre- 
sponding profiles (to visualize them, we used a Gaussian 
profile with amplitude u n at site n). In addition, we 
found two non-symmetric solutions 



0(7, A) 



-7A 2 ± ,/7M 4 -4A 2 (l + A 2 ) 2 
2A 2 (1 + A 2 ) 
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Without loss of generality, we restrict to the case \a\ < 1 . 
The sign "+" in (Q} corresponds to the asymmetric so- 
lution (|a oav | 7^ 1) that exists only for 7 > 4, bifur- 
cating from the antisymmetric mode [see Fig. HJa) at 
A ~ 8]. Once this asymmetric solution appears, the anti- 
symmetric one becomes unstable [a standard linear sta- 
bility analysis [22j was performed and full (dashed) lines 
indicate stable (unstable) solutions in Fig. [Tj. Fig. [TJa) 
shows a monotonic increment of the asymmetric power 
up to some maximum value. All this branch (full red 
line) is stable. However, after achieving this maximum, 
the branch changes its curvature (power decreases) and 
the solution gets unstable until it fuses with the anti- 
symmetric branch [this unstable branch corresponds to 



the sign "— " in Eq. Q]. For saturable systems, this 
always unstable non-symmetric solution is called inter- 
mediate solution 0-@. The case where two fundamen- 
tal solutions are simultaneously stable, sharing the same 
Hamiltonian value, was initially suggested [B| as a van- 
ishing Peierls-Nabarro barrier [2J, [24[. However, recent 
works HH have shown that there is a nonzero effec- 
tive energy barrier, which also considers the IS. Surpris- 
ingly, the dimer model also shows this phenomenology 
which is fundamental to understand deeply the proper- 
ties of nonlinear saturable arrays. Examples for some 
profiles are sketched in Fig. [TJb). The intermediate and 
asymmetric modes correspond both to non-symmetric so- 
lutions of model JT]); they have quite similar profiles and 
would be only identified by directly observing their un- 
stable/stable dynamical propagation. 



A. Effective potential 

Now, we go deeper into the dynamical properties of this 
model by computing an effective potential 0"H, HE HH ■ 
Model (dJ is integrable, therefore the effective potential 
can be obtained analytically. First of all, by fixing the 
Power P, we define the center of mass as x = u 2 /P (x = 



or 1 



P 



and x = 0.5 



P/2). 



Then, as we arc considering a defocusing nonlincarity, 
we study staggered solutions (a < 0) and express the 
amplitudes as 



mi = ±\J P(l — x) and u 2 = ^fVxP . 

With these expressions inserted into the Hamiltonian ([3]) , 
we get 
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Figure 2: (Color Online) Effective potential for 7 = 10 and 
P = 0.15 (dotted), P = 15 (dashed), P = 19 (full) and P = 26 
(full thick) . Filled circles correspond to stationary solutions. 



The critical points dH/dx = represent different sta- 
tionary solutions. A first one corresponds to x ant = 0.5 
(q = —1), a solution existing for all level of powers P [see 
blue line in Fig.QJa)]. There are four additional solutions 



-± — a/P 2 -2 7 2 
2 2P V 



4(P+l)±2 7v / 7 2 -4(P- 



1), 
(5) 

where asymmetric modes correspond to the sign "+" in 
front to the inner square root. They are symmetrically 
located - to the right and to left - from the antisym- 
metric solution (x an t = 1/2). Asymmetric solutions ex- 
ist only for 7 > 4, bifurcating from the antisymmetric 
mode at power P m j„ = 7 — 2 — (7 — 2) 2 — 4. Unsta- 
ble IS correspond to the sign "— " in front to the inner 
square root in ((5|) . They exist above the power threshold 
7 — 2 + y/ ("f — 2) 2 — 4. The asymmetric and the 
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intermediate solutions exist up to an upper power value 
given by P up = ( 7 /2) 2 - 1. (For 7 = 10, P min = 0.254, 
P th = 15.746, and P up = 24 [see Fig. Ufa)]). IS also ex- 
ist for 7 > 4; i.e, the intrinsic saturable phenomenology 
would be only observed above some critical nonlinearity. 

Fig. [2] shows the effective potential H(x,P, 7) versus 
the center of mass, for 7 = 10 and for four different 
level of power (curves have been normalized for compar- 
ison). Below the P m in (dotted line), the potential looks 
as a typical potential well with a minimum located at 
Xant- In the range {P m in,Pth} the effective potential is 
cubic- like (dashed line); i.e, the antisymmetric solution 
is unstable (maximum) while the asymmetric solution is 
stable (minimum). The saturable nature of this model 
manifests for powers above Pth (full line). The antisym- 
metric and the asymmetric solutions become, both, a lo- 
cal minimum and, therefore, simultaneously stable. The 
unstable symmetry-broken IS appears as a maximum in 
this potential, located in between the stable solutions. 
Therefore, at Pth the IS has a center of mass xjs — x an t 
that evolves in the direction of an asymmetric configu- 
ration when approaching P up . Finally, above P up , both 
non-symmetric solutions disappear, and the antisymmet- 
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Figure 3: (Color Online) Density Plot of u x ut" versus \a\ 
and P for 7 = 10. Full and dashed white lines correspond to 
the asymmetric and intermediate solutions, respectively. 



ric mode is the only critical point in the effective potential 
(full thick line). This is a direct consequence of the sat- 
urable nature of the system, where the nonlinear term 
vanishes for powers above a critical value (P > P up ). 



B. Numerical propagation 

Now, in order to test these stationary properties 
and their dynamical consequences, we numerically study 
model ([T|) by considering the general initial condition 
lti(0) = A and u 2 (0) = aA, with -1 sC a < 0. This 
input condition allow us to excite all different staggered 
solutions. Fig.[3]shows the output center of mass, defined 
as Xout = \u2{z m ax)\ 2 / P , measured after a given propa- 
gation distance z = z max (similar density maps are ob- 
tained for different propagation distances). Purple-blue 
colors (Xout ^ 0.2) represent solutions localized close to 
site n = 1; green colors (0.4 < x out < 0.6) represent an- 
tisymmetric profiles; orange-red colors (x ou t > 0.8) rep- 
resent solutions localized close to site n = 2. For P < 1, 
asymmetric profiles (\a\ ~ 0) do not correspond to any 
stationary solution and the light just oscillates between 
the two sites (see the appearance of multiple colors as 
an indication of strong oscillation, i.e. switching [15|). 
On the other hand, for \a\ ~ 1 we see a greener color 
that indicates a small oscillation in the vicinity of the 
antisymmetric solution {x ou t = 0.5), the only stationary 
solution at this power regime. At the bifurcation point 
power (when H changes its shape from dotted to dashed 
in Fig. [J) , both solutions are quite similar and the anti- 
symmetric mode is slightly unstable. Then, in the range 
P G {1, 15} for \a\ < 0.4, the light is well trapped at the 
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vicinity of site n = 1 . This is an indication of the exci- 
tation of an asymmetric stationary state (sec full white 
line in Fig. |3]). For \a\ — > 1, there is an oscillation of the 
light in the interval x ou t € {0,0.5}, as expected for an 
unstable antisymmetric configuration (see dashed line in 
Fig. [2J. In the region of power ps {16, 24}, a very inter- 
esting behavior is observed: we found a simultaneously 
stable dynamical propagation of asymmetric (|cv| « 0.15) 
and antisymmetric (|a| ~ 1) input profiles. The constant 
value of x (constant color) when increasing P, indicates 
a stable dynamical evolution. This constitutes a dynam- 
ical and indirect evidence of the existence of the inter- 
mediate solutions as a fundamental entity for saturable 
systems: if both staggered solutions are stable simultane- 
ously, an extra unstable intermediate solution must ex- 
ist. The dashed white line in Fig. [3] corresponds to the 
Intermediate solution and shows a clear connection in 
parameter space between the asymmetric and the anti- 
symmetric staggered solutions. For P > 24, we observe 
an oscillation of energy between sites 1 and 2 for the 
input condition < |a| < 0.8, while an antisymmetric 
configuration is stable for \a\ > 0.8. All the dynamical 
results are in perfect agreement with the stationary pic- 
ture sketched in Figs. Q] and [U including the region of 
multi-stability and the existence of the IS. 



IV. EXPERIMENTAL RESULTS 




PM 



BS 



L2 
AM 




M 



cw 532nm 



CCD 



MO 



SNC 



MO 



BS 



Figure 4: (Color Online) Experimental setup for observing 
the multistable regime in a saturable nonlinear coupler. 

To verify our theoretical and numerical predictions we 
use the experimental setup sketched in Fig. @] A cw laser 
with a wavelength of 532 nm propagates through a phase 
mask (PM) covering half of the beam along the trans- 
verse direction x (coinciding with the crystallographic 
c-axis of the sample). So the phase relation of the left 
and right half of the beam can be switched to be ei- 
ther in-phase or out-of-phase. With a 4f imaging sys- 
tem, composed by lenses LI and L2, the beam is imaged 
onto a double-hole amplitude mask (AM). By using a 
microscope objective, the beam is injected into a sat- 
urable nonlinear coupler (SNC) fabricated by titanium 



in-diffusion on an a;-cut lithium niobate substrate doped 
with iron. The end facet of the sample is monitored 
by a high-resolution CCD camera. The length of our 
sample along the propagation z-direction is 18 mm (~ 5 
diffraction lengths) where the waveguide channels are 4.0 
/zm wide with a separation of 2.2 /xm. Our photovoltaic 
samples have a nonlinearity which grows exponentially 
in time, ^(t) = 7(1 — exp[— t/r}), where r is the dielec- 
tric response time [27j ■ In order to reach a steady-state, 
saturation of the photovoltaic nonlinearity is required, 
which typically occurs, in our samples, for t ^ 25 min- 
utes. Fig. [5] shows some examples of the evolution of the 
center of mass in time. We see how some input condi- 
tions propagates in a very stable way for the whole mea- 
surement period (black-full, orange-full and black-dashed 
lines) while the orange-dashed curve tend to stabilize for 
t > 25. 
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Figure 5: (Color Online) Center of mass evolution at the 
sample output for asymmetric (full lines) and antisymmet- 
ric (dashed lines) input conditions, for input powers lOOnW 
(black) and 500nW (orange). The filled rectangular area de- 
notes the region where we experimentally observe (define) an 
antisymmetric configuration. 

We repeat the experiment for several input powers and 
for the asymmetric and the antisymmetric input profiles, 
and average the center of mass values for the last 10 min- 
utes of each experiment. Compiled results are shown in 
Fig. [5] with a fairly good agreement between the theoret- 
ically predicted phenomenology and what is observed in 
direct experiments. The initial conditions are indicated 
by dots connected with dashed lines. The deviation of the 
output values (dots connected by full lines) from the ini- 
tial conditions can be understood as the degree of stabil- 
ity. In addition, the error bars (obtained in the averaging 
process) also indicates how stable is the profile: a smaller 
bar indicates a dynamically stable profile while a larger 
bar means a stronger oscillation around some minimum 
(stationary solution). For low level of power, we observe 
that both solutions are essentially stable. This coincides 
with our analysis in the region of the first bifurcation 
point, where the appearance of the asymmetric solution 
weakly destabilizes the antisymmetric mode. Then, for 
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Figure 6: (Color Online) Averaged center of mass at the sam- 
ple output vs input power for a LiNb03 saturable coupler. 
Red and blue symbols, connected by full lines, correspond to 
the experimentally measured value of x ut for an asymmet- 
ric and antisymmetric input condition (symbols connected 
by dashed lines), respectively. The lower row shows some 
(non-averaged) experimental output profiles for the indicated 
powers with antisymmetric (ant) and asymmetric (asy) input 
profiles. (The filled area denotes the same that in Fig. [5}. 

powers in the region ~ {300, 1300}nW, a typical cubic- 
like picture is observed: the asymmetric solution is sta- 
ble while the antisymmetric one is not. However, above 
~ 1300nW a saturable phenomenology emerges. The sat- 
urable nature of the nonlinearity allows oscillation of the 
stability properties, with a richer dynamics when com- 
paring with usual cubic systems. For this sample, the 
region ~ {1300, 2300}nW corresponds to a regime where 
the two fundamental nonlinear solutions of this problem 
become simultaneously stable (as an example see the pro- 
files for 1500nW). Therefore, we experimentally observe 
a multi-stable regime as an indirect evidence of the ex- 
istence of intermediate solutions in a saturable dimer. 
This observation gives a strong support to the theory 
(and model) developed for this type of lattices. Above 
some given power (~ 2300nW), the only stable solution 



is the antisymmetric one, as shown in Fig. [6] The large 
bar for the last asymmetric (red) point indicates that the 
system has saturated and thadit the effective potential 
has a shape like the thick line in Fig. [5] It is impor- 
tant to mention that our experimental output profiles in 
Fig- EJ- low-row have a staggered phase, what is evident by 
observing the zero amplitude in between the two waveg- 
uides. This gives us an extra support for relating the 
observation of the multi-stable regime with the existence 
of an unstable IS, because the observed profiles are well 
connected in phase space. 

V. CONCLUSIONS 

In conclusion, we have observed, for the first time, a 
multi-stable regime of fundamental modes in a nonlinear 
saturable coupler. A complete map of nonlinear solutions 
has been constructed, including their stability properties, 
and effective potential. Numerically, we have determined 
the regions where the fundamental profiles are expected 
to be stable, showing an excellent agreement with the 
developed theory. We fabricated a nonlinear saturable 
coupler in LiNbC»3 and observed stable propagation of 
fundamental modes for intermediate level of power. This 
constitutes the first experimental evidence, in any phys- 
ical system, of the existence of intermediate solutions 
in discrete nonlinear lattices. Moreover, for larger pow- 
ers we observed the absence of the asymmetric solution 
as a clear indication of the saturation of the nonlinear- 
ity. All these results support strongly recent theoretical 
and numerical developments based on DNLS-like mod- 
els in several contexts of physics and open, as a direct 
consequence, new opportunities of research on nonlin- 
ear discrete systems. In addition, the understanding of 
small systems becomes very crucial when thinking on the 
implementation of photonic lattices for realistic applica- 
tions. 
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